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When a quasi-particle Andreev reflects from a normal- 
superconducting (N-S) interface, the phase of the outgo- 
ing excitation is shifted by the phase of the supercon- 
ducting order parameter [1]. Consequently if a phase 
coherent normal conductor is in contact with two super- 
conductors with order parameter phases 0i, 4)2, trans- 
port properties will be oscillatory functions of the phase 
difference 4> = (pi — <p2- Spivak and Khmel'nitskii [2] 
argued that the ensemble average of the electrical con- 
ductance of a diffusive structure would have the form 
{G{4>)) = A + Bcos{2(j)) were A and B are positive con- 
stants of magnitude less than the quantum of conduc- 
tance 2e^/h. Thus the ensemble averaged conductance 
is predicted to have a meLximum at zero phase and to 
have a fundamental periodicity of tt. Later it was noted 
[3] that the conductance G of a single sample should have 
a periodicity of 2tt, with an amplitude < 2e^//i, although 
no comment on the possibility of a zero phase minimum 
was offered. The theories of references [2,3] were formu- 
lated before exact formulae [4,5] for the electrical conduc- 
tance of N-S-N structures were available and therefore, 
to avoid significant supercurrents arising from Andreev 
scattering, were restricted to the regime where the dis- 
tance between the superconductors is greater than the 
thermal coherence length Lt = yhD/kBT . 

The aim of this Letter is to address two crucial ques- 
tions raised by recent experiments on such interferom- 
eters [6-10]. The first concerns the nature of the zero 
phase extremum in G{(f>), which in the theory of [2] 
and all experiments using normal contacts is found to 
be a maximum, whereas in the low temperature limit 
M < Lt, theory [11] predicts that minima are allowed. 
A second concerns the amplitude of oscillation, which in 
the experiments [7,8,10] is found to be smaller than or 
of order 2e^//i, but in those of references [6,9] is sev- 
eral multiples of 2e^/h. If the latter are not an arti- 
fact of the experiments, then superconductivity enhanced 
quasi-particle interference effects (SEQUINs) of this kind 
should be present in the low temperature theory of refer- 
ences [4,5,11]. In what follows we demonstrate that large 
amplitude SEQUINs are indeed obtainable from exact so- 
lutions of the Bogoliubov - de Gennes equation and high- 
light conditions under which giant oscillations should be 
observable. Remarkably we predict that the amplitude 
vanishes for both very dirty and very clean samples and 
that to obtain a large effect, a degree of normal scattering 
must be introduced, in order that an approximate sum 
rule is broken. Further we predict that the experiments 
carried out to date are far from optimal and that am- 
plitudes of oscillation many orders of magnitude greater 
than e'^/h are attainable. 

For simplicity we consider the zero temperature limit, 
where the electrical conductance between two normal 



reservoirs can be written [4,5] (in units of 2e^/ft,), 

2{RaR'a - TaT'^) 



G = n + Ta 



Ra + R'a+Ta+ 
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In this expression, Ro,To {Ra,Ta) are the coefficients 
for normal (Andreev) reflection and transmission for 
zero energy quasi-particles from reservoir 1 , while i?g , Tq 
(i?^, T^) are corresponding coefficients for quasi-particles 
from reservoir 2. If each of the external leads connecting 
the reservoirs to the scattercr contains N open channels, 
these satisfy Ro + Tq + Ra+Ta = R'a + U + K + T'a=N 
and Tq -\- Ta = Tq + T^. Furthermore, in the absence 
of a magnetic field, all reflection coefficients are even 
functions of 0, while the transmission coefficients sat- 
isfy T^(0) = To{-4), T^i4) = Ta{-4>). Consequently 
on quite general grounds, in the absence of a field, G is 
predicted to be an even function of 4 [12]. 

Figure 1 shows three examples of interferometers, for 
which results are presented below. Each has two super- 
conducting regions with definite phases (pi and (j)2, in 
contact with a normal region (shown shaded). In each 
case the scattering region is connected to normal, exter- 
nal current carrying leads, with N conducting channels. 
In figures la and lb, a normal barrier (shown black) is 
placed at the N-S interface. In what follows we show 
numerical results obtained from a tight binding model, 
with diagonal elements and nearest neighbour hopping 
elements of magnitude unity. In the external leads e; = 
and in the single line of sites forming the barrier = £(,. 
In a disordered region of the sample, is chosen to be 
a random number uniformly distributed over the inter- 
val ±.W . As discussed in reference [5], the conductance 
is obtained by first evaluating the quantum mechanical 
scattering matrix and then evaluating the zero tempera- 
ture conductance formula (1). The transfer matrix codes 
used [5] are extremely versatile and can be used to ana- 
lyze arbitrary geometries, with multiple contacts. 

For the structure of figure la, figure 2 shows numeri- 
cal results for giant oscillations in the electrical conduc- 
tance G{(f)). In the absence of a barrier (e{, = 0), the 
amplitude of oscillation (in units of 2e^/h) is negligible 
compared with unity, whereas for £{, = 1, 2, 3 a large am- 
plitude oscillation is present. As the barrier strength 
increases, the amplitude of oscillation initially increases 
to a value of order Ne^ /h, before decreasing in propor- 
tion to the zero phase conductance. These results show 
that at intermediate barrier strengths, the relative am- 
plitude as well as the absolute amplitude is optimised. 
Figure 3 shows results for the phase periodic conduc- 
tance of structures lb (sohd hue) and Ic (dashed line). 
For the structure lb in figure 3 we have presented results 
for the most favourable barrier strength, in the absence 
of disorder. Introducing normal disorder or changing the 
barrier strength decreases the amplitude of oscillation. 
In the case of structure Ic in figure 3 there is no barrier. 
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but a disorder comparable to the experiments of refer- 
ence [9] has been used. We have examined the structures 
in figures lb and Ic for a variety of disorders and barrier 
strengths and in no case have found an ampUtude which 
is more than a few percent of N2e'^ /h. 

The crucial role of normal scattering in optimising this 
effect can be understood through a multiple scattering 
description of a N-S interface, which captures the essen- 
tial physics of interferometers. Consider first an idealisa- 
tion of the structure of figure la, in which the distance M 
between the superconductors vanishes and therefore for 
a long enough sample there is no quasi-particle transmis- 
sion. In this limit the total resistance reduces to a sum of 
two measureable boundary resistances and in what fol- 
lows, we therefore focus attention on the left-hand bound- 
ary conductance [13] 

N 

Gb{cP) = 2Ra = 2Tr r„rt = (^a)y (2), 

where {Ra)ij ~ l(^a)u P is the Andrecv reflection proba- 
biltity from channel j to channel i. As in equation (1), 
the Andreev reflection coefficient is of the form Ra = 

-^diag ^-^ofl — diag 

where i?diag — lZi=li^a)ii and i?off-diag 

is the remaining contribution from inter-channel scatter- 
ing, i?off-diag = J2i^j=li^a)ij- 

In the absence of disorder, for M = and </> = 0, trans- 
lational symmetry in the direction transverse to the cur- 
rent flow guarantees that i?off-diag = 0. For the system of 
figure la, with no barrier, no disorder and M = 45, figure 
4b, shows the behaviour of the coefficients {Ra )ij for i j 
and demonstrates that even for finite M, off-diagonal 
scattering at ^ = is negligible. This figure leads us 
to a second observation, namely that even for non-zero 
(f), almost all of the off-diagonal coefficients arc negligibly 
small and that a given channel i couples strongly to at 
most one other channel j. Consequently in the absence of 
disorder, a multiple scattering description involving pairs 
of channels captures the essential physics. 

Consider a normal barrier to the left of a N-S interface. 
Particles (holes) impinging on the normal scatterer are 
described by a scattering matrix Spp, (shh), and those 
arriving at the N-S interface by a reflection matrix p, 
where 

g —I '''pp ^pp\ p — ( I 

V hv ^'ppj ' \ Php Phh J ' 

The elements of s and p are themselves matrices de- 
scribing scattering between open channels of the external 
leads. For an ideal interface, where Andreev's approxi- 
mation is valid [1] , ppp and phh are negligible and in what 
follows will be set to zero. As a consequence, php and 
Pph are unitary and one obtains [14] = t'f^f^phpM~ptpp, 
with Mpp = 1 — r'ppPphr'f^i^php- In contrast with the anal- 
ysis of [14], where php is proportional to the unit ma- 
trix, the interference effect of interest here is contained 



in the fact that p^p induces off-diagonal scattering. Sub- 
stituting Ta into equation (2) and taking advantage of 
particle- hole symmetry at E = 0, yields 

G = 2Tr(TQ-ir(gt)-i) (3) 

where Q = p*^ -|- {r')ppPph{r')lp, with T = tppt^p the 
transmission matrix of the normal scattering region. This 
multiple scattering formula for the boundary conduc- 
tance is valid in the presence of an arbitrary number of 
channels and in any dimension. Notice that if T is equal 
to the unit matrix, Q = p^f^ and therefore G = 2N, irre- 
spective of the phase periodic nature of pph- This demon- 
strates that at a clean interface, whatever the phase, 
the approximate unitarity of pph yields the sum rule 
^diag+.Roff-diag = N and therefore the conductance is in- 
dependent of 4>. More generally, whenever normal reflec- 
tion {Ro) and Andreev transmission (Tq) are negligibly 
small, unitarity imposes the sum rule Tq + Ra = N and 
since in this limit equation (1) reduces to G = Tq + Ra, 
the amplitude of oscillation must vanish. 

Equation (3) is very general and makes no assumption 
about the nature of matrices pph and Spp. We now intro- 
duce a two-channel model in which pph is chosen to be an 
arbitrary two dimensional unitary matrix. In the absence 
of disorder, tpp and rpp are diagonal and therefore the 
only interchannel coupling is provided by pph ■ Substitut- 
ing these matrices into equation (3), yields an expression 
for fa involving a single phase 9, whose value is a linear 
combination of phase shifts due to normal reflection at 
the barrier. Andreev reflection at the N-S interface and 
the phase accumulated by an excitation travelling from 
the barrier to the interface. In what follows, we present 
results for the average value of (Ra)^j, obtained by in- 
tegrating over 6. For a given value of (/), once the nor- 
mal barrier transmission coefficients Ti and T2 of the two 
channels are chosen, {Ra)ij is completely determined. 

For the structure of figure la, figures 4a and 4b show 
numerical results for the diagonal {Ra)ii and off-diagonal 
coefficients {Ra)ij {i j) respectively. Notice that 
at zero phase, most of the diagonal coefficients {Ra)ii 
are close to unity, although a small number of order 
N\A\/Ep, where |A| is the order parameter magnitude 
and Ep the Fermi energy, are suppressed, due to a break- 
down of Andreev's approximation for low angle scattering 
[15]. This slight breakdown of Andreev's approximation 
yields a small amplitude oscillation even in the absence 
of normal potential scattering, but as emphasized by fig- 
ure 2, the fractional amplitude is negligible. Figures 5a 
and 5b show corresponding results for the diagonal and 
off-diagonal coefficients in the presence of a barrier. At 
(f> = 0, there is no coupling between the channels and the 
scattering properties are those of N independent chan- 
nels, each with a barrier transmission coefficient Tj. The 
spectrum of the coeffcients dcpcmds in detail on the shape 
of the barrier. The inset of figure 5a shows the bound- 
ary conductance G{4>) obtained by summing the curves 
in figures 5a nad 5b. Figures 5c and 5d show results from 
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the two channel evaluation of equation (3), obtained by 
choosing 10 pairs of transmission coelRcients Ti, T2, with 
T2 = O.2T1. Clearly the qualitative features of the ex- 
act simulation are reproduced by this simple two channel 
analysis. 

As shown in figure 2, in the absence of zero phase inter- 
channel scattering, the zero phase extremum is a mini- 
mum. Prom figure 4, it is clear that the nature of the 
extremum is the result of a competition between diago- 
nal Andreev reflection coefficients, which exhibit a zero 
phase maximum and off-diagonal coeflacients which pos- 
sess a zero phase minimum. Analytically wc find that, 
in the absence of disorder, the second derivative of two 
channel conductance satisfies cPG/dcj)^ > 0, for all barrier 
strengths. In contrast, figure 6 shows numerical results 
for the structure of figure la, with M = 45, M' = 50, 
M" = 15, but with the barrier replaced by a disordered 
normal square of width 30 sites. This shows that re- 
placing the barrier by a disordered region, causes the 
extremum of the off-diagonal coefficients to switch from 
a minimum to a maximum. In addition, the channels 
no longer couple in pairs and a complete multi-channel 
scattering description is needed. It is noted [3,5] that 
changing the microscopic impurity configuration shifts 
G{(j)) by an amount of order e'^/h, so in the presence of 
giant oscillations, the nature of the zero phase extremum 
is insensitive to such changes. 

In summary, wc have demonstrated that giant SE- 
QUINs are obtainable from exact solutions of the Bo- 
goliubov - de Gennes equation and can be observed by 
breaking a crucial sum rule. Remarkably, the structure 
of figure Ic used in the experiments of [9] shows only a 
small percentage effect, which nevertheless, in view of the 
large number of channels in these experiments, yielded an 
amplitude greater than 2e'^/h. We predict that a more 
optimal structure is obtained in the presence of a normal 
barrier at the interface and that in metallic samples, with 
very large A^, the amplitude of oscillation could be orders 
of magnitude larger than 2e^//i. 

Since carrying out this work, we became aware of a re- 
lated paper on diffusive conductors [16], where it is pre- 
dicted that {Ra)ii ^ {Ra)ij {i 7^ j) ■ This phenomenon 
of giant diagonal scattering (i.e. backscattering ) was of- 
fered as the origin of large oscillations in G((^). In this 
Letter, we predict (sec eg. fig. 4) that giant oscillations 
can occur even when diagonal and off-diagonal scattering 
coeffiecients are comparable. 
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FIG. 1. Three possible interferometers, each with two su- 
perconducting regions of width M" . In figures la and lb, the 
superconductors are separated by a distance M and in Ic by a 
distance 3M. In figure la, the scattering region is connected 
to normal, external current carrying leads, of width M+2M", 
in figures lb and Ic of width M . In figures la and lb, a nor- 
mal barrier (shown black) is placed at the N-S interface. The 
current flows from left to right between external reservoirs 
with potentials /ii and /i2. In the tight binding model used 
in the numerical simulations, the barrier comprises a line of 
sites with diagonal elements €» = £(,. 

FIG. 2. Numerical results for the conductance G of the 
structure of figure la, with M = 45, M' = 50, M" = 15, 
and number of open channels N = 45. Results are shown 
for barrier potentials tt, = 0, 1, 2, 3. The number adjacent a 
given curve is the correponding value of ej. 



FIG. 3. Figure 3 shows results for the conductance G 
of the structure of figure lb (solid line: left scale) with 
M = M' = M" = 10, W = Q, and Ic (dashed line: right 
scale) with M = M' = M" = 10, W = 0. A = 10 in both 
cases. For the structure lb results are shown for a barrier 
potential €6 = 2. 
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FIG. 4. Figure 4a shows numerical results for the diagonal 
Andreev reflection coefficients {Ra)ii of the structure of figure 
la, with M = 45, M' = 50, M" = 15, iV = 45, and no barrier 
present, (ej, = 0.) Figure 4b shows corresponding results for 
the off-diagonal coefficients {Ra)ij with i =^ j. 

FIG. 5. Figure 5a shows numerical results for the diagonal 
Andreev reflection coefficients {Ra)u of the structure of figure 
la, with M = 45,M' = 50, M" = 15, N = 45, and barrier 
potential, = 2. Figure 5b shows corresponding results for 
the ofF-diagonal coefficients {Ra)ij with i ^ j. Figures 5c and 
5d show analytical results from a two channel calculation. 

FIG. 6. Numerical results obtained from a tight bind- 
ing model of the structure of figure la, but with the barrier 
replaced by a disordered region of length 30 sites. In those 
simulations, M = 45, M' = 15, M" = 50, N = 45, Ao = 0.1, 
and the disorder is W = 2.8. 
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